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Abstract 

We report on the results of a numerical simulation concerning the low-lying 
spectrum of four-dimensional J\f = 1 SU(2) Supersymmetric Yang-Mills (SYM) 
theory on the lattice with light dynamical gluinos. In the gauge sector the tree-level 
Symanzik improved gauge action is used, while we use the Wilson formulation in the 
fermion sector with stout smearing of the gauge links in the Wilson-Dirac operator. 
The ensembles of gauge configurations were produced with the Two-Step Polynomial 
Hybrid Monte Carlo (TS-PHMC) updating algorithm. We performed simulations on 
large lattices up to a size of 24^ • 48 at /3 = 1.6. Using QCD units with the Sommer 
scale being set to tq = 0.5 fm, the lattice spacing is about a ~ 0.09 fm, and the 
spatial extent of the lattice corresponds to 2.1 fm. At the lightest simulated gluino 
mass the spin- 1/2 gluino-glue bound state appeared to be considerably heavier than 
its expected super-partner, the pseudoscalar bound state. Whether supermultiplets 
are formed remains to be studied in upcoming simulations. 
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1 Introduction 



In recent years supersymmetric theories have aroused increasing interest in elementary 
particle physics. The supersymmetric extension of the Standard Model with M = 1 
supercharge is considered to be an interesting candidate for a quantum field theory with 
phenomenological relevance in the near future. Supersymmetry (SUSY) is an essential 
ingredient also for other models beyond the Standard Model. 

The M = 1 Supersymmetric Yang-Mills (SYM) theory is the minimal supersymmetric 
extension of the SU(Ai'c) gauge theory describing self-interactions of gauge fields A^, cor- 
responding to the gluons [g). The supersymmetric partners of the gluons are described by 
spin-1/2 Majorana fermion fields A" (a = 1, . . . , — 1), the gluinos (g). Compatibility 
of SUSY with gauge invariance requires that the gluinos transform in the adjoint repre- 
sentation of the gauge group. This theory describes the interactions between gluons and 
gluinos. The Lagrangian of Euclidean SYM theory in the continuum, including a SUSY 
breaking mass term, reads 

1 1 — 7T? ~ — 

^SYM = + -x'^i.iv.xr + ^ A-A"^ , (1) 

where denotes the gauge covariant derivative in the adjoint representation. The gluino 
mass term introduces a soft breaking of supersymmetry. 

In the low-energy regime the interactions become strong. Arguments based on the 
low-energy effective Lagrangian approach [2] predict the occurrence of non-perturbative 
dynamics like confinement and spontaneous chiral symmetry breaking in SUSY gauge 
theories. Confinement is realised by colourless bound states. Since both gluons and 
gluinos transform according to the adjoint representation, bound states can be built by 
any number of at least two gluons and gluinos. In the case where the last term in Eq. ([T]) 
is switched off {rrig = 0), an anomalous global chiral symmetry U(1)a is present. This 
symmetry is equivalent to the i?-symmetry in supersymmetric models. The anomaly does 
not break the global chiral symmetry completely and a discrete subgroup Z2Nc remains. 
As in the case of QCD, the discrete chiral symmetry is expected to be spontaneously 
broken to Z2 by the non- vanishing value of the gluino condensate (AA). The consequence 
of this spontaneous breaking is the existence of N^. degenerate ground states with different 
orientations of the gluino condensate. 

SYM is also equivalent to QCD with a single quark flavour {Nf = 1 QCD) in the 
limit of a large number of colours {N^ — )■ 00), where the Majorana spinor is replaced by a 
single Dirac spinor in the antisymmetric representation of the gauge group [3j. The latter 
model is also object of investigation by our collaboration [4]. 

Since confinement occurs in low-energy SYM, standard analytical methods like per- 
turbation theory fail and non-perturbative methods are required. This motivates the 
introduction of the lattice formulation of SYM. The first lattice formulation of SYM suit- 
able for numerical simulations has been proposed by Curci and Veneziano [5]. It it based 
on the Wilson discretisation, which proved to be successful in lattice QCD computations 
in spite of its known limitations. First non-perturbative investigations of SYM on the 
lattice using this formulation have been performed by [B] in the quenched approximation, 
and by the DESY-Miinster-Roma collaboration with dynamical fermions; see Ref. [7| for 
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a review, and references [HI El [IDl [HI IE] • SUSY is broken explicitly by the lattice discreti- 
sation. Additionally, in the Wilson approach the mass term and the Wilson-term break 
both chirality and SUSY explicitly. Both symmetries are expected to be recovered in the 
continuum limit by tuning the relevant bare mass term to its critical value corresponding 
to a massless gluino (mg = 0), and the gauge coupling towards zero. 

In recent years, simulations of A/" = 1 SYM on the lattice using Ginsparg- Wilson 
fermions with good chiral properties, such as domain wall fermions, have been initiated 
[131 [m [151 US] • For large lattice volumes and small lattice spacings these formulations 
require, however, a significantly larger amount of computing resources than the Wilson 
formulation. The gain of no need for tuning the position of the zero gluino mass point 
does not compensate by far the advantage of Wilson fermions. 

In the past, investigations of the gluino dynamics have been performed using the 
Two-Step Multi Bosonic (TSMB) algorithm [T7]. This algorithm was developed in the 
framework of the DESY-Miinster collaboration. Recently, the Two-Step Polynomial Hy- 
brid Monte Carlo (TS-PHMC) algorithm [TH] has been developed and implemented for 
SYM. This algorithm offers more efficiency and improvements compared to the TSMB 
algorithm and allows us to collect higher statistics and to simulate small gluino masses 
rrig in this study. Furthermore, due to available computer resources we simulated the 
theory on volumes with extension larger than 2 fm, which is expected to be the minimally 
required volume for spectroscopic studies. 

The main purpose of this work is to continue the project of the DESY-Miinster col- 
laboration for the simulation of A/" = 1 SU(2) SYM. We present new accurate results 
obtained with the newly used TS-PHMC algorithm and improved actions. 

The most important characteristics of the theory is the mass spectrum of bound states, 
for which the low-energy effective theories predict a reorganisation of the masses in two 
massive Wess-Zumino supermultiplets at the SUSY point [HIS], where the soft breaking 
vanishes. The introduction of a small gluino mass removes the mass degeneracy between 
the supermultiplet members. In the lower supermultiplet the ordering of the states with 
increasing mass is: scalar glueball 0"*"+, spin- 1/2 gluino-glueball (xl), pseudoscalar glue- 
ball 0""*". The ordering is reversed in the higher supermultiplet which contains: adjoint 
pseudoscalar meson a-rj', spin-1/2 gluino-glueball (xh), adjoint scalar meson a-fo- 

The plan of this paper is as follows: in the next Section we review the lattice formula- 
tion and describe the simulation details. Section [3] is devoted to the static quark potential 
and the determination of the scale. Methods for the determination of the masses of bound 
states are described in Section (H In Section \5\ the results on the spectrum are collected 
and discussed. Finally, we conclude our findings in Section [HI 

2 Lattice formulation oi J\f = 1 SYM theory 

The Curci-Veneziano action of A/" = 1 SU(2) SYM theory on a lattice, S = Sg + Sg, 
contains the usual plaquette gauge field action Sg, and a fermionic action Sg for the 
gluino. The gauge action Sg can be extended to a more general form which includes, 
besides the usual (1 x 1) Wilson loop plaquette term, (1 x 2) Wilson loops of perimeter 
six. We employ the tree-level improved Symanzik (tlSym) gauge action, given for SU(Ai'c) 
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colour group by 

with the normahsation condition cq = 1 — 8ci. The bare gauge couphng is related to 
the lattice parameter /3 by the usual relation /3 = 2Nc/gQ. For the tlSym action we have 
ci = -1/12 [19]. 

The gluinos are represented by Majorana fermions A" in the adjoint representation. 
They satisfy the Majorana condition 

X = X^ = CV, (3) 

where C = 7072 is the charge conjugation matrix in the spinorial representation. 

In the gluino sector, the Wilson formulation for fermions proposed in [5J introduces 
the Wilson term proportional to r, which is an irrelevant term in the continuum limit. 
We set the Wilson parameter to r = 1. The fermion part Sg of the action is then given 
by 

^^ = ^EM^)A(x)-|EE[M^+/i)V;(x)(l+7,)A(x)+A(x)Vf(x)(l-7,)A(x+/i^^ (4) 

X X fl 

where k is the bare hopping parameter which encodes the bare gluino mass k = (2mg + 
8)~^. The real orthogonal matrices V^{x) are the gauge links in the adjoint representation: 

[V,{x)f ^ 2Tr[Ul{x)T^U,{x)T''] = [V;{x)f = [V-\x)r , (5) 

where are the generators of SU(Ai'c) satisfying 2Tr(T'^T^) = 5"^. In case of SU(2) one 
has T°- = ^a"' with the Pauli matrices a". 

The links Ux^fi in the fermion action can be replaced by stowi-smeared links [20j. 
This has the advantage that short range topological defects of the gauge field and the 
corresponding small eigenvalues of the fermion matrix are removed. Both the tlSym 
gauge action and the stout smeared links in the fermionic part of the lattice action are 
introduced in order to accelerate the approach to the continuum limit as /3 — )■ cxd. 

The stout smeared links are defined by 

Ui]l ^ exp { ^ (fi.,, - Ql^) - ^ Tr (^],,, - fit } . (6) 

Here ?7^^^ denotes the original "thin" gauge links, and 

nx,,^pUl^C.,, (7) 

with the sum of "staples" 

Cx,fi = E/ {yx+ji,v^x+v,yUx,v + Ux-ii-{-[i,uUx-ii,^iUx-C',u) ■ (^) 



4 



p is an arbitrary parameter which we fix in this work to p = 0.15. In principle, the 
smearing defined by the above equations can be iterated several times, but then the 
fermion action becomes extended over a larger region on the lattice. We prefer to keep 
the action well localised and hence only perform a single smearing step. 
Writing the gluino action as 

•^9 ~ 9 ^ '^^ \Qyx^x , (9) 

Q is the non-hermitian fermion matrix or lattice Wilson-Dirac operator for Dirac fermions 
in the adjoint representation. Using relation ([H]), the fermion action can be rewritten in 
terms of the antisymmetric matrix M = CQ. Integration of the fermionic variables yields 
the Pfaffian of M, 

ypAe"^« = Pf(M), (10) 
whose absolute value equals the square root of the fermion determinant: 

|Pf (M)| = y^det(M) = y^det(Q) . (11) 

Effectively, this corresponds to a fiavour number Nf = 1/2. In the Wilson setup, det(Q) 
and det(M) are always real and positive, but the Pfaffian Pf (M) can become negative 
even for positive gluino masses. 

In our numerical simulations we include the dynamics of the gluino by the Two- 
Step Polynomial Hybrid Monte Carlo (TS-PHMC) |18j algorithm with fiavour number 
Nf = 1/2. This has the consequence that only the absolute value of the Pfaffian is taken 
into account in the updating of the gauge field configuration. The sign of the Pfaffian 
has to be included in a reweighting step when calculating expectation values. It can be 
shown that the sign of the Pfaffian is equal to the sign of the product of half of the 
doubly degenerate negative real eigenvalues of Q. For positive gluino masses sufficiently 
far away from zero, a negative sign of the Pfaffian rarely occurs in the updating sequence 
and therefore in this situation a sign problem does not show up. Approaching the limit 
of vanishing gluino mass we monitor the sign of the Pfaffian and take it into account by 
reweighting. It turned out that only in our runs D and Dg (see Tabled] below) a noticeable 
number of configurations with negative sign occured; the highest fraction being in point 
Dg, where they amount to 3 % of all configurations. The effect of the negative signs on 
the particle masses turned out to be negligible. 

The parameters of the J\f = 1 SYM on the lattice are the lattice gauge coupling /3 
and the fermionic hopping parameter k. Similarly to QCD, the mass term proportional 
to rrigfi breaks chirality explicitly. In the present case it also breaks the supersymmetry. 
A massless gluino, rrig = 0, is obtained by tuning the bare mass term to its critical value 
(^g,o — ^ TTT^c) or equivalcutly k — Kc- 

In order to study questions related to supersymmetry, one has to approach the critical 
value of the hopping parameter k = Kc corresponding to zero gluino mass. This tuning 
problem can be solved rather easily by means of the adjoint pion mass m^r. This is 
the pion mass in the corresponding theory with two Majorana fermions in the adjoint 
representation. It is obtained from the exponential decay of the connected part of the 
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Table 1: Algorithmic parameters for TS-PHMC runs with tlSym gauge action at /3 = 1.6. 
Runs labelled with subscript s have been performed with Stout-links. Nconf is the number of 
configurations produced, Tq is the Sommer scale parameter, am^^ is the adjoint pion mass 
in lattice units, and Mj. is the dimensionless quantity Mr = (romTr)^ used to estimate the 
gluino mass. In Mr the values of r^/a extrapolated to k,c have been used. 



Run 




(3 


K 


Nconf 


ro/a 




Mr 


A 


16^ 


32 


1.6 


0.1800 


2500 


2.9(1) 


1.3087(12) 


45.6(4.2) 


B 


16^ 


32 


1.6 


0.1900 


2700 


3.3(1) 


1.0071(12) 


27.0(2.5) 


C 


16^ 


32 


1.6 


0.2000 


10847 


4.242(87) 


0.5008(13) 


6.68(62) 


D 


16^ 


32 


1.6 


0.2020 


6947 


5.04(26) 


0.221(12) 


1.30(19) 


A 




48 


1.6 


0.1980 


1480 


3.885(63) 


0.6415(13) 


11.0(1.0) 


B 




48 


1.6 


0.1990 


1400 


4.16(12) 


0.5759(17) 


8.83(82) 


C 


243 


48 


1.6 


0.2000 


6465 


4.33(19) 


0.4947(13) 


6.52(61) 


As 


24^ 


48 


1.6 


0.1500 


370 




0.9469(38) 


28.69(89) 


Bs 


243 


48 


1.6 


0.1550 


1730 


4.324(39) 


0.5788(16) 


10.72(33) 




243 


48 


1.6 


0.1570 


2110 


5.165(88) 


0.3264(23) 


3.41(11) 


Ds 


243 


48 


1.6 


0.1575 


2260 


5.561(99) 


0.2015(93) 


1.30(13) 



pseudoscalar meson propagator, see below. The pion is not a physical particle in the 
spectrum of the SYM theory, but it can be unambiguously defined in a partially quenched 
framework. To determine the pion mass is rather easy, in fact it is the easiest mass to 
determine. As will be detailed in Sec. 14.41 the behaviour of the pion mass-squared is 
very closely linear as a function of 1/k in the entire range of gluino masses of interest. 
On the basis of arguments involving the OZ I- approximation of SYM [I], the adjoint pion 
mass is expected to vanish for a massless gluino. Therefore it is enough to perform two 
simulations on relatively small lattices at relatively large gluino masses, from which Kc 
can be obtained by a linear extrapolation. Proceeding to larger lattices and smaller gluino 
masses, this estimate can be continuously improved without any further simulations. In 
this way the interesting range of hopping parameters k < for the investigation of the 
particle spectrum can be determined. 

The values of the gauge coupling parameter /3 can be fixed by investigating the static 
potential of an external charge in the fundamental representation and extracting the 
Sommer scale parameter ro/a [2T], as discussed in Sec. El In analogy with QCD, we set 
the value of tq by definition to tq = 0.5 fm. In this way we can use familiar QCD units 
for physical dimensionful quantities. 

As a measure for the gluino mass we define the dimensionless quantity Mr = (rom,r)^; 
which is expected to be proportional to the gluino mass. 

A summary of the simulation parameters is given in Table [TJ The simulations are 
performed on 16'^ ■ 32 and 24^ ■ 48 lattices. Extrapolated to Kc the lattice spacing amounts 
to a ~ 0.097 fm for the unstout ensembles and a ~ 0.088 fm for the stout ones, see below. 
The lattice extension L ~ 2.1 — 2.3 fm is expected to be large enough to allow control 
over finite volume effects on the masses of the bound states. 
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An issue in lattice simulation is the lightness of the dynamical fermions which leads 
to slowing down of the update algorithms. The TS-PHMC algorithm turned out to be 
very efficient in producing short autocorrelations among the gauge configurations. For 
instance, in the stout-smeared runs on a 24^ ■ 48 lattice the integrated autocorrelation 
of the average plaquette (which belongs to the worst quantities from the point of view 
of autocorrelations) did always satisfy rf^^'^ < 10. The lightest adjoint pion mass in our 
simulations was about 440 MeV. Simulations for smaller gluino masses and/or finer lattice 
spacings are going on presently. 



3 Static potential and physical scale 

Analogy with QCD suggests that the colour charge is confined in SYM, so that the 
particle states are colour-singlets. Moreover, SYM is expected to confine static quarks 
as in pure Yang-Mills theory: the static quark-antiquark potential can not be screened 
by the dynamical gluinos transforming in the adjoint representation, and a non- vanishing 
string tension arises at large distances. 

The numerical results for the static potential V{r) for runs A-D are shown in Fig. [TJ 
The linear behaviour at large quark-antiquark separations is compatible with a non- 
vanishing string tension. 




Figure 1: The static quark potential in A/" = 1 SU(2) SYM. The solid 
lines are fits to the data. 

From the behaviour of the static potential at intermediate distances it is possible to 
determine the lattice scale, a well-known procedure in lattice QCD. The scale can be 



7 



characterised by the Sommer parameter ro [21] defined by the relation 



2 ^ 
dT 



1.65. (12) 



ro 



In the SYM model the string tension could in principle also be used to fix the scale, but 
the Sommer scale parameter is more convenient from the numerical point of view. The 
numerical results for vq/q are reported in the sixth column of Tabled! Following the 
analogous procedure in QCD in a mass independent renormalisation scheme, we extrap- 
olate these data to vanishing adjoint pion mass. For the runs with thin links, where we 
combine data from the two volumes, we obtain ro/a = 5.16(24), and for runs with stout 
links ro/a = 5.657(85). Using ro = 0.5 fm this corresponds to a = 0.097 fm (thin links) 
and a = 0.088 fm (stout links), respectively. The physical size of the simulated boxes is 
therefore in these units L ~ 1.5 — 2.3 fm. 



4 Spectrum of low-lying bound states 

For the investigation of the spectrum of low-lying bound states we concentrate on the 
operators employed for the construction of the low-energy Lagrangians of [Tj and [2]. 
These are expected to dominate the dynamics of SYM at low energies. Previous experience 
on the determination of low-lying masses is reported in |8j and [TU] . We investigate spin-0 
gluino-gluino bilinear operators (adjoint mesons), a spin-1/2 mixed gluino-glue operator 
and spin-0 glueball operators. In some cases smearing techniques such as APE [22] and 
Jacobi smearing [23] are applied in order to increase the overlap of the lattice operator 
with the low-lying bound state. 

4.1 Adjoint mesons 

The adjoint mesons are colourless states with spin-parity 0^ and 0^, composed of two 
gluinos. In analogy to fiavour singlet states in QCD we denote the former a-rj' and the 
latter a-/o, where the prefix a indicates "adjoint". The associated projecting operators 
are the gluino bilinear operators Omeson = where F = 1 or P = 75, respectively. The 
resulting propagator consists of connected and disconnected contributions: 

C{xo - yo) = Cconn(a;o - Vo) - C'disc(a;o - Vo) 

= yEm^Qxl^Qyi]) - ^E(Tr[rg-L]Tr[PQj), (13) 

^ X ^ X 

where (• • • ) denotes the average over the gauge sample and Vg = L^. 

The exponential decay of the connected part defines the adjoint pion mass 771^^- This 
quantity, even if not associated to a physical state of SYM, can be used to determine 
the gluino mass, as mentioned in Sec. [21 Indeed, according to arguments involving the 
OZI-approximation of SYM [T], the adjoint pion mass is expected to vanish for a massless 
gluino and the behaviour ml oc rrig can be assumed for light gluinos [H [TU] . 
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As is well known in simulations of QCD, the numerical evaluation of the disconnected 
propagators is rather demanding. We employ here two alternative methods, the Stochastic 
Estimators Technique (SET) in the spin dilution variant |24j, and the Improved Volume 
Source Technique (IVST) [22]. As in QCD, the disconnected diagrams are intrinsically 
noisier than the connected ones and dominate the level of noise in the total correlator. 



Lattice: 16 X 32 k=0.2 
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8 10 



12 14 



16 



Figure 2: Connected and disconnected pieces and the total time-slice 
correlation function of the adjoint pseudoscalar a-rj'. 

In the pseudoscalar channel a reasonable signal-to-noise ratio is obtained, allowing the 
extraction of the mass from the mass fit on most samples. As an example, in Fig. [2] we show 
the result for the a-rj' correlator for run C together with the two different contributions. 
Examples for the effective masses as a function of the time separation t are shown in 
Fig. El 

In the scalar channel the extraction of the mass is complicated by the presence of a 
vacuum expectation value for the projecting operator ~ (AA). This allowed a relatively 
precise determination of the a-/o mass only for the samples with stout smearing, which 
give a better signal. Two examples for the effective mass in this channel are shown in 
Fig. m For the future we plan the application of variance reduction techniques for a more 
precise computation of the disconnected diagrams. 

4.2 Scalar glueball 

As for the adjoint mesons, we investigated the scalar glueball masses also in both parity 
channels. In order to improve the signal we applied in this case APE smearing with 
the variational method [26j. For the positive parity glueball 0"*" we adopted the simplest 
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Figure 3: Effective mass of the pseudoscalar a-t]'. The horizontal hne 
represents the result from a one- mass-fit. 



a-L Lattice: 24^ 48 
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Figure 4: Effective mass of the scalar a-/o. The horizontal line repre- 
sents the result from a one-mass- fit. 

interpolating operator built from space-like plaquettes: 

Ogiue,+(a;) = Ttc[Ui2{x) + ^723(0;) + Uu{x)] . 
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For the negative parity state we considered the eight-hnk operator proposed in |H] . How- 
ever, the signal obtained in this case was too poor to obtain an estimate of the mass. 
Therefore we restricted the analysis to the positive parity channel in the following. 

Also here, as for the scalar a-/o, the gauge samples generated with stout links generally 
turn out to give better results for the glueball masses. In Fig. \5\ two examples of the 
effective masses are reported together with the results from one-mass fits with minimal 
time-distance tj. 
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Figure 5: Effective mass and results from one- mass fits for the scalar 
glueball 0+. 



4.3 Gluino-glueballs 

The gluino-glueballs (gg) are spin 1/2 colour singlet states of a gluon and a gluino. They 
are supposed to complete the Wess-Zumino supermultiplet of the adjoint mesons [I]. For 
this state we adopt the lattice version of the gluino-glue operator TrcfFcrA] P, where the 
field-strength tensor Ff^^{x) is replaced by the clover-plaquette operator Pfj,u{x) P [TU]: 

(^f.(^) = E4'Tr.[P.,(x)A^(a:)]. (15) 

Here only spatial indices are taken into account in order to avoid links in the time- 
direction. The clover-plaquette operator, having the correct behaviour under discrete 
parity and time reversal transformations, is defined as 

p,Ax) = ^j:{ui:^{x) - u;^:>Hx)) (16) 
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with 



U^'Jix) = Ul{x)Ul{x + u)U,{x + (i)U,{x) (17) 

U^K^) = Ul{x)U,{x -v + fi)U,{x- i))Ul{x -v) (18) 

UfX^) = - f>)Uu{x -V- fl)Ulix - - f,)Ulix - A) (19) 

Ul'Jix) = U,ix - jl)Ul{x - jl)Ul{x + 0- fl)U,ix) . (20) 

The full correlator of the gluino-glue operator, 

Cgixo - yo) = E E «Tr[f/,,(x)a1Q;i, ^,,,Tr[f/,,(l/)aVf;0 , (21) 

S,y i,j,k,l 

is a matrix in Dirac space with two independent components 

Cg^ixo - yo) = C.ixo - yo)5"^ + C,,ixo - yohf , (22) 

with Ci = TTD[Cgg]/4: 0,11(1 C2 = Tr£i[7oCgg(x)]/4. We see agreement in the masses 
extracted from each component, see Fig. [HI For the final estimates we choose the time 
antisymmetric component Ci, which appears to provide better plateaus. We apply APE 
smearing for the links and Jacobi smearing for the fermion fields in order to optimise the 
signal-to-noise ratio and to obtain an earlier plateau in the effective mass. 
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Figure 6: Effective mass of the gluino-glueball gg. 
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4.4 Massless gluino limit 

Of high interest in lattice simulations of SYM is the point corresponding to a massless 
gluino, where supersymmetry is expected to emerge in the continuum limit. With Wilson 
fermions this point must be located by a tuning procedure due to the additive renormal- 
isation of the bare gluino mass. 

The subtracted gluino mass can be determined in different ways. It can be directly 
obtained from the study of lattice SUSY Ward-Identities (Wis) as discussed in [H]. We 
have implemented the determination of the necessary operators for the Wis. Apart from 
confirming the smallness of lattice corrections to the Wis, consistent with 0{a) effects, 
they give results for the gluino mass up to a renormalisation factor. On the other hand, 
the point of vanishing gluino mass can be estimated in an indirect way from the vanishing 
of the adjoint pion mass. Indeed, as mentioned above, the pion mass squared (am^)^ is 
expected to vanish linearly with the (renormalised) gluino mass. 

Both the Wis and adjoint pion mass methods give consistent estimates of the critical 
hopping parameter Kc corresponding to vanishing gluino mass. As an example, in Fig. [7] 
we show the the gluino mass and the pion mass squared as a function of 1/k for the 
runs with unstout links on the 24^ ■ 48 lattice. Both clearly show a linear behaviour. 
The linear extrapolations to vanishing gluino mass give = 0.2027(4) from the Ward 
identities and k^^^ = 0.20300(5) from the pion mass. Similarly, for the runs with stout 
links on the 24'^ ■ 48 lattice we obtain = 0.15883(85) from the Ward identities and 
K^^^ = 0.15793(4) from the pion mass, which agree within errors. 




4,95 5 
1/K 



Figure 7: The gluino mass from the SUSY Ward identities (left panel) 
and the pion mass squared (right panel) as function of the inverse hop- 
ping parameter The critical value Kc is indicated by the asterisk 
symbol. 
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5 Spectrum of bound states 



The masses of the hghtest bound states of low-energy A/" = 1 SYM determined in this 
work are collected in Table [2] and a graphic representation is shown in Fig. [SI 

Table 2: Results for the low-lying bound state masses of Af = 1 SU(2) SYM from the 
various runs. The masses are given in lattice units. 



Run L K, 


a-T]' a-fo gg glub. 0++ 


A 16 0.1900 
B 16 0.1800 
C 16 0.2000 
D 16 0.2000 


1.3115(67) 2.229(80) 1.862(21) 1.291(13) 
1.0396(72) 1.27(18) 1.546(14) 1.156(51) 
0.5425(71) 0.931(53) 0.982(10) 0.941(18) 
0.361(60) 0.87(10) 0.7532(96) 0.819(19) 


A 24 0.1980 
B 24 0.1990 
C 24 0.2000 


0.675(18) 1.15(12) 1.1456(82) 
0.6215(86) 1.314(32) 1.0789(95) 
0.536(24) 0.863(81) 0.9895(70) 0.781(24) 


As 24 0.1500 
Bs 24 0.1550 
Cs 24 0.1570 
Ds 24 0.1575 


1.0114(82) 1.07(16) 1.302(14) 
0.614(23) 0.964(70) 0.9559(48) 1.079(92) 
0.416(29) 0.467(93) 0.7250(56) 0.582(61) 
0.327(30) 0.351(85) 0.682(30) 0.389(90) 



The masses in Fig. [H] are multiplied by the extrapolated value of the Sommer scale 
parameter and plotted as a function of the squared adjoint pion mass for {rorriT^)'^ < 12. 
The lightest simulated adjoint pion mass is about 440 MeV in our units. The vertical line 
in Fig. [H] indicates the massless gluino limit where SUSY restoration is expected up to 
0{a) effects. The physical extent of the lattice is 1.5 — 2.3 fm. 

The bound state masses appear to be characterised by a linear dependence on (romTr)^, 
in accordance with the prediction of [2Z!- An extrapolation of our data with stout links 
(points Bs to Ds), which have better numerical quality than the unstout ones, to the 
massless gluino limit yields the numbers in Table [31 

Table 3: Bound state masses in physical units extrapolated to the massless gluino limit. 





a-T]' a-fo 


99 


glub. 0++ 


m [MeV] 


670(63) 571(181) 


1386(39) 


721(165) 



The gluino-glueball {gg) with a mass of about 1386 MeV turns out to be considerably 
heavier than the a-rj' with a mass of 670 MeV. Furthermore, the masses of the scalar glue- 
ball and the scalar meson a-fo are near the mass of the pseudoscalar a-rj'. The behaviour 
of scalars is compatible with mixing between 0"*" glueball and a-fo- The pattern of scalar 
masses suggests a lower supermultiplet, while the spin- 1/2 candidate remains heavier up 
to the smallest simulated gluino mass in this simulation, and also after extrapolation to 
Kc- Whether this outcome is a discretisation artefact or a physical effect, as claimed in 
[25] . should become clear in future studies at finer lattice spacings. As the data at small 
gluino mass are preliminary, it would be premature to make judgements about this point. 
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Spectrum of N=1 SU(2) Super- Yang-Mills theory on the lattice 
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Figure 8: Low-lying bound state masses of A/" = 1 SU(2) SYM as a 
function of the adjoint pion mass squared in physical units. The blue 
symbols represent the extrapolations to the massless gluino limit. 



6 Conclusion 

In this work first quantitative results on the low-energy spectrum of A/" = 1 supersymmetric 
Yang-Mills theory are obtained. Physical volumes larger than 2 fm have been simulated, 
which is the volume usually required for spectrum studies in lattice gauge theory. The 
comparison of masses on different volumes in otherwise same conditions reveals negligible 
finite size effects at least for moderate gluino masses. We have collected higher statistics 
and have used efficient dynamical algorithms such as TS-PHMC, which is suitable for light 
fermion masses. In addition, the supersymmetric Ward identities and other observables 
like the confinement potential have been investigated. 

From the results of the mass spectrum the question of the gluino-gluino and gluino- 
glueball mass splitting remains open. It can only be answered by further simulations 
allowing an extrapolation to the continuum limit. 
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